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ABSTRACT 

We calculate the flux of cosmic gamma rays expected from the annihilation of neutralinos in the Galactic 
halo. Our calculation of the annihilation cross section to two photons improves the existing calculations by 
inclusion of exact one-loop diagrams for the amplitudes involving Higgs boson and chargino states as well as 
those involving fermion and sfermion states. A survey of supersymmetric parameter space shows that numerous 
models would be observable at the 3cr level with an air Cerenkov telescope with an exposure of 10 4 m 2 yr. 
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1. INTRODUCTION 



For quite some time, it has been known that only about a tenth of the mass 
of most galaxies, including our own, is luminous, and that the rest is composed 
of some sort of dark matter [1]. The nature of this nonluminous material is 
unknown, although there are quite convincing arguments that it must be non- 
baryonic. One of the most promising of many candidates for the dark matter 
is the neutralino [2] [3] , a linear combination of the supersymmetric partners of 
the photon, Z°, and Higgs bosons. It has been suggested by numerous authors 
that if neutralinos populate the Galactic halo, then monoenergetic gamma rays 
produced by neutralino annihilation in the halo could provide a plausible avenue 
toward discovery of such dark-matter particles [4] . In this paper we re-examine 
this proposal. 

The goal of this work is to provide results for the cross section for anni- 
hilation of neutralinos to two photons which include all of the contributions 
at one loop to the amplitude for a neutralino in any given minimal supersym- 
metric extension of the standard model. We generalize previous calculations of 
the amplitude for annihilation through quark-squark loops [5] [6] [7] and Higgs- 
chargino loops [5] to arbitrary neutralino and squark masses and compositions. 
We also include the recent approximate calculation of Bergstrom and Kaplan [8] 
of annihilation through TT /± -boson-chargino loops and improve it by including 
subleading logarithmic terms. 

We then use these expressions in a survey of supersymmetric parameter 
space to assess the possibility of discovering dark-matter neutralinos in the 
Galaxy via observation of cosmic gamma rays produced by neutralino annihi- 
lation. 

In the following, we will present our result for the cross section, discussing 
the importance of the various contributions. Then we will give an estimate of 
the signal rates which are implied. 
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FIG. 1. Feynman diagrams for neutralino annihilation to photons. 

2. CROSS SECTION 

The Feynman diagrams for annihilation to two photons are shown in Fig. 1. 
The diagrams fall into three categories. The similarity between the diagrams of 
classes a) and b) indicates that the corresponding amplitudes can be written in 
terms of the same basic functions arising from the loop integrations. Further- 
more, these functions are precisely those which appear in the calculation of the 
cross section for neutralino annihilation to two gluons presented in Ref . [9] . The 
third class of diagrams are those with W bosons in the loop, and those ghost 
diagrams which are related to them; this gauge invariant set of diagrams has 
been discussed in Ref. [8]. By a choice of non-linear gauge, the calculation was 
reduced significantly. In the limit m x ± >m x o (where m x ± and m x o are masses 
of the chargino x ± and neutralino x° i respectively) , which is always appropri- 
ate when considering the neutralino as the lightest supersymmetric particle, the 
amplitude reduces to a single three-point integral, as given in Ref. [8]. 
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Neutralinos in the halo move with velocities negligible compared with the 
speed of light, so annihilation occurs in the s wave only. Therefore, the ampli- 
tudes depend only on the outgoing photon momenta and polarizations, and the 
amplitudes can be written in the form 



A = ■^e(k 1 ,k 2 ,e 1 ,e 2 )A, 



(2.1) 



where fc, and e, are the momenta and polarizations of the outgoing photon pair. 
The total amplitude will be a sum of three parts to be discussed in turn below, 
A = Ajj + Arc + Aw- Given this amplitude, the cross section is 



a 2 m 2 ~ „ 



(2.2) 



First, consider the amplitudes related to the two gluon amplitudes, i.e., 
those coming from fermion-sfermion loop diagrams shown in class a) in Fig. 1. 
Define the following functions, arising from the loop integrals: 
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(2.4) 



where tua and mz are the masses of the pseudoscalar-Higgs and Z bosons, 
QAxx an d 9z X x are tne couplings of neutralinos to the pseudoscalar-Higgs and 
Z bosons [9], and 
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Note that this expression for F(a, b, S, D) corrects a typographical error in 
Ref. [9]. 

The amplitude in terms of these functions is 

+ T— f(^,^,s ~,d -)\, 

f j j J 

, 2 (2-6) 

Im^^-.E^K-^ln^j-T^,^,^) 

+ m 2 L> 7 1 m 2 - ' m 2 - ' [ ' 

where the sum on / is over quarks and leptons, and the sum on / is over the 
squarks and sleptons. Here, gAff are the pseudoscalar-Higgs-fermion couplings 
[9], /3f = (1— m^/m^) 1 / 2 , Qf is the electric charge of /, c/ is a color factor which 
equals 3 for quarks and 1 for leptons, and 6{x) is the Heaviside step function. 
The function I(x) is given in Eq. (2.14) of Ref. [9]. The couplings of the 
fermions and sfermions are defined by = A 2 - + B 2 ~ and D ^ = A 2 - — B 2 ~, 
with 

A ff = \( X ffo + W ffo)> {27) 

B ff = \ ( v ;-/u - w ff-o) ■ 

The fundamental couplings X' and W are the couplings of left-handed and 
right-handed fermions, respectively, to sfermions and neutralinos, as defined in 
Refs. [9] and [10]. 

Next consider the amplitude from diagrams involving intermediate charged- 
Higgs bosons and charginos, i.e., those in class b). We find 



{9 9 
X X 




(2.8) 



where the sum is over the two charginos. Here S m = (g/ cos 9 w)(0'^ m — 0^ m ) 
is the coupling of the Z to a chargino pair, and h Ax ± x ± = <?(sin (3 Q mm + 
cos (3S mm ) is the coupling of the A boson to a chargino pair. The couplings of 
the charginos to neutralinos are defined by 



s m = \[(Q'L) 2 + (QZ) 2 ] 



o' L n' R 



(2.9) 



Om- 



Definitions and more details of the couplings O' and Q 1 can be found in Refs. 
[2] and [10]. 

Next consider the amplitude from diagrams with intermediate W bosons, 
i.e., those in class c). In Ref. [8], the imaginary part of this amplitude was 
calculated exactly, 



ImAw = 7T 9{m 2 - m w )C m (3 W In ] + ^ w 
*— J 1 — p 



w 



(2.10) 



where (3w = (1 — m^/m^) 1 / 2 . The coupling is given by 



(2.11) 



where 0§ m and 0^ m are given in Refs. [2] and [10]. The real part of the am- 
plitude was evaluated in a leading- logarithm limit in Ref. [8], using dispersion 
relations, giving a result proportional to In 2 (a). As pointed out there, it is an 
excellent approximation to contract the chargino propagator to a point, cor- 
responding to the limit m x ±>m x . Evaluating the real part of the subsequent 
three-point amplitude, we find 
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(2.12) 



(2.13) 



and n is the renormalization point; note that n appears because of the contrac- 
tion of the chargino propagator. We choose [i = mw, specifying the running 
coupling at the electroweak scale; this gives a subleading logarithmic contri- 
bution to the amplitude. It is very useful to have an asymptotic expansion 
for this result, extending the leading-logarithm expansion of Ref. [8] into the 
subleading terms. We find 



r> / \ { -±\r? a +1 + \ a 



4r a -> 0, 

(2.14) 

t-> / \ f — ^lna— 1 a — > 0, 
#2 (a) ~ < 1 i 

1-2 "Si a ^°°- 

The leading logarithmic behaviour of this amplitude agrees with that obtained 
in Ref. [8] . It should be pointed out that these expressions are not strictly valid 
for m x < mw- However, since we are interested in high energy gamma rays, 
this case is not of interest to us; we henceforth assume that m x > mw- 

Now we can consider the behaviours of each of these terms. The size of 
the fermion/sfermion loop contribution is sensitive to the gaugino fraction of 
the neutralino and to the masses of the sfermions; a sufficiently Higgsino-like 
composition of the neutralino will suppress these contributions. As pointed out 
by Bergstrom and Kaplan [8] , annihilation through loops containing W bosons 
[class c)] may be significant. 

The Higgs-chargino loop contribution can be important for several reasons. 
First note that the couplings which appear in these diagrams are all naturally 
of the order of the gauge couplings, so that there is no special suppression, 
as can sometimes occur in couplings of neutralinos to fermions and sfermions. 
Furthermore, the masses of the intermediate particles are comparable, and this 
matching of scales provides some enhancement of the loop integral. In fact, 
it is important to note that the amplitude as written contains a pole at the 
point b = 1 + a, or m 2 + m 2 H ± = m 2 ± . This divergence occurs because we 
have ignored the widths of the intermediate particles. In actuality these widths 
are large; for chargino masses greater than 100 GeV the width of the chargino 
is approximately T x ± ~ 0.1m x ±. So the pole is spurious, and the amplitude 



must be mollified in the region b ~ 1 + a. Note that previous treatments of this 
contribution [5] will not be reliable in the common case that the masses of the 
particles are comparable. 

The TV-loop contribution is almost always important in the case that the 
neutralino is heavier than the W and primarily Higgsino [8]. This is because 
this contribution depends most strongly on the lightest-chargino mass and this 
is usually not too much larger than the neutralino mass. 



3. SIGNAL RATES 

We follow Urban et al. [11] and consider the signal from pointed obser- 
vation of the Galactic center with an atmospheric Cerenkov telescope (ACT). 
This is potentially a promising method for observation of high-energy gamma 
rays from neutralino annihilation. 

The flux of gamma rays from neutralino annihilation, from a window of 
solid angle Aft aimed at the Galactic center, may be written [4], 

^ * 2 x lO" 11 cm" 2 sec" 1 (p^) 2 f h J^I^^^l {An /sr) ^ 

(3.1) 

where /3° 4 is the local halo density in units of 0.4 GeV cm~ 3 , and 0.3 < /halo is 2 
is a parameter which reflects uncertainty in modeling the galactic halo. Eq. 
(3.1) is obtained assuming an isothermal halo with a density profile suitable 
for accounting for the observed rotation curve. If the cosmological neutralino 
density is too small to account for the halo dark matter, neutralinos should still 
gather in galactic halos, and the annihilation cross section in this case should 
generally be larger than that in the case where neutralinos are the dark matter. 
Therefore, there may be an observable gamma-ray signature of neutralinos even 
if they exist but are too few to account for all the halo dark matter. To account 
for this, we take the halo density (for ft x h 2 < 1 as required by the age-of-the- 
Universe constraint) to be fl x h 2 /0.25, where fl x is the cosmological neutralino 
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FIG. 2. Minimum exposure required for a 3<r detection of gamma rays from neu- 
tralino annihilation in the Galactic center, versus mass of the neutralino, for the 
survey of supersymmetric parameter space discussed in the text. The symbols indi- 
cate which of the three amplitude contributions dominates the cross section; triangles 
indicate that the fermion-sfermion diagrams dominate, diamonds indicate that the W 
diagrams dominate, and circles indicate that the Higgs-chargino diagrams dominate. 



abundance, and h is the Hubble parameter in units of 100 km sec -1 Mpc -1 . 
We note that our results will not depend sensitively on this prescription. 

The natural width of the two-photon peak is small, and the background 
below the peak is controlled by the instrumental resolution. The most impor- 
tant background for energies below 1 TeV comes from misidentified charged 
particles [11]; the background from diffuse cosmic gamma rays [6] contributes 
only slightly to the total background in this regime. Following Ref. [11], we 
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find a background flux of 

An 



AE 

b b = (1.2 x lO^secTV^): " 



10" 3 sr 



L100 GeVJ 



(3.2) 



100 GeV 

where AE is the energy resolution on the peak, and p ~ 3.3 is the spectral index 
for the background, in this case dominated by misidentified electrons [11]. 

In Fig. 2 we plot the exposure required for a 3<r detection of gamma rays 
from neutralino annihilation in the Galactic center for a variety of models. The 
parameter ranges which generated these models were taken to be 100 GeV < 
M 2 < 800 GeV, 200 GeV < \x < 800 GeV, 2 < tan/3 < 20, 300 GeV < m A < 
600 GeV, and 200 GeV < m,- < 800 GeV. Here M 2 and /i are the gaugino 
mass parameters, tua is the mass of the pseudoscalar Higgs particle, and m 9 - is 
a common squark mass parameter (which differs from the actual squark masses 
due to mixing terms, which were included). The grand unification condition 
was assumed. Models were cut from the plot if they violated known bounds 
from e + e~ physics, if they gave Higgs masses in volation of current limits, or if 
they were inconsistent as models for neutralino dark matter (for example, we 
obviously require that the lightest supersymmetric particle is a neutralino). 

Triangles, circles, and diamonds indicate models where the fermion- 
sfermion, Higgs-chargino, and TU-boson diagrams, respectively, dominate the 
amplitude. According to Fig. 2, the fermion-sfermion diagrams are most often 
important, but there are indeed regions of parameter space where the Higgs- 
chargino and W loops are dominant. These results indicate that numerous 
supersymmetric models could be probed by an atmospheric Cerenkov detector 
with an area £>(10 4 m 2 ). 



4. ACKNOWLEDGMENTS 

G.J. thanks J. Buckley for conversations about high-energy gamma rays. 
M.K. acknowledges the hospitality of the Theory Group at CERN where part of 
this work was completed. M.K. was supported at the I.A.S. by the W. M. Keck 
Foundation and at Columbia University by the U.S. Department of Energy 
under contract DE-FG02-92ER40699. G.J. was supported by the U.S. D.O.E. 
under contract DE-FG02-85ER40231. 



9 



References 



[1] For recent reviews of dark matter and its detection, see V. Trimble, Ann. 
Rev. Astron. Astrophys. 25, 425 (1989); J. R. Primack, B. Sadoulet, and 
D. Seckel, Ann. Rev. Nucl. Part. Sci. B38, 751 (1988); Dark Matter in the 
Universe, eds. J. Kormendy and G. Knapp (Reidel, Dordrecht, 1989). 

[2] H. E. Haber and G. L. Kane, Phys. Rep. 117, 75 (1985). 

[3] J. Ellis, J. S. Hagelin, D. V. Nanopoulos, K. A. Olive, and M. Sred- 
nicki, IVud. Phys. B238, 453 (1984); K. Griest, M. Kamionkowski, and 
M. S. Turner, Phys. Rev. D 41, 3565 (1990); M. Drees and M. M. Nojiri, 
Phys. Rev. D 47, 376 (1993); K. A. Olive and M. Srednicki, Phys. Lett. 
B230, 78 (1989); K. A. Olive and M. Srednicki, Mid. Phys. B355, 208 

(1991) ; J. McDonald, K. A. Olive, and M. Srednicki, Phys. Lett. B283, 80 

(1992) ; K. Griest, Phys. Rev. D 38, 2357 (1988); FERMILAB-Pub-89/139- 
A (E); Phys. Rev. Lett. 61, 666 (1988). 

[4] M. Kamionkowski, to appear in The Gamma-Ray Sky with Compton GRO 
and SIGMA, proceedings of the Winter School, Les Houches, France, Jan 
5-Feb 4, 1994. 

[5] G. F. Giudice and K. Griest, Phys. Rev. D 40, 2549 (1989). 
[6] S. Rudaz, Phys. Rev. D 39, 3549 (1989). 
[7] L. Bergstrom, Phys. Lett. B225, 372 (1989). 

[8] L. Bergstrom and J. Kaplan, preprint USITP-94-03, hep-ph/9403239 
(1994). 

[9] M. Drees, G. Jungman, M. Kamionkowski, and M. Nojiri, Phys. Rev. D 
49, 636 (1994). 

[10] K. Griest, G. Jungman, and M. Kamionkowski, to be submitted to Physics 
Reports. 

[11] M. Urban et al., Phys. Lett. B293, 149 (1992). 



10 



